Entropic commensurate-incommensurate transition 
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The equilibrium properties of a minimal tiling model are investigated. The model has extensive 
ground state entropy, with each ground state having a quasiperiodic sequence of rows. It is found 
that the transition from the ground state to the high temperature disordered phase proceeds through 
a sequence of periodic arrangements of rows, in analogy with the commensurate-incommensurate 
transition. We show that the effective free energy of the model resembles the Frenkel-Kontorova 
Hamiltonian, but with temperature playing the role of the strength of the substrate potential, and 
with the competing lengths not explicitly present in the basic interactions. 

PACS numbers: PACS numbers: 64.60.De , 61.44.Fw, 61.44.Br 



Tilings provide a simple means to model systems with 
both simple and complex ground states. As statistical 
mechanical models, they include such systems as the 
Ising and Potts models, as well as other models with 
discrete spin variables [ I] . The class of tilings, however, 
is much larger|9j, and includes the various non-periodic 
Wang tilings 0, the quasiperiodic Penrose tilings, the 
asymptotically isotropic Pinwheel tiling [10], and the re- 
cently discovered generalizations of the Rudin-Shapiro se- 
quence, which are neither periodic nor quasiperiodic |2"0] . 

The statistical mechanical behavior of tiling models is 
rich, and, as yet, largely uncharted. Apart from results 
which may be transcribed from discrete spin models, only 
a very small number of systems have been studied. First, 
a model based on the Amman set of 16 Wang tiles with 
quasiperiodic ground states was studied by several au- 
thors [111 IT31 IT?] . It appears that this model undergoes 
a continuous phase transition from a disordered state to 
a quasiperiodic phase, and its non-equilibrium behavior 
was studied in the context of spin-glasses. A variation of 
the model allowing more complicated interactions and 
vacancies shows a first order transition |4J. Hierarchi- 
cal tilings [14] have also been studied, and a very recent 
model[5] possesses limit-periodic ground states which un- 
dergo a series of phase transitions where motifs of ever 
larger scales order as the temperature is lowered. Finally, 
we note some recent studies |181 HO] on models with large 
number of degenerate disordered ground states aimed at 
studying glasses. 

In this paper we study the equilibrium behavior of a 
model based on the 13-tile Kari-Culik (KC) set |7J [JO] 
of Wang tiles, both numerically and analytically. The 
KC set is the smallest known aperiodic set - it is the 
smallest set of tiles which can tile the plane, but not pe- 
riodically. Allowed juxtapositions of tiles are enforced 
by matching rules, and these in turn induce a Hamil- 
tonian: Every matching rule violation is penalized by a 
positive energy, while allowed matchings have zero en- 
ergy. In what follows, we shall denote the energy cost of 
mismatching adjacent vertical or horizontal edges by J x 
or J y , respectively. 

We shall argue that the equilibrium behavior of this 
system is analogous to the Frenkel-Kontorova (FK) 
model|6 of the commensurate-incommensurate (CI) 
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Figure 1: The 13 KC tiles. The seven upper tiles are type A, 
and the six lower tiles are type B. Note that |, and 0' are 
considered different markings. 



transition. The FK model describes a chain of masses 
which are connected by springs and subjected to a pe- 
riodic substrate potential. It exhibits rich behavior due 
to competition between these two interactions, each of 
which favors ordering with a different wavelength. How- 
ever, in marked contrast with the FK model, where the 
favored length scales are present in the Hamiltonian, in 
the KC model they emerge spontaneously, from only 
nearest-neighbor interactions. Moreover, the role of sub- 
strate potential strength in the FK model is played by 
the temperature in the KC model. In this sense, the KC 
system exhibits an entropic CI transition. 

As indicated in Figure [TJ the 13 KC tiles may be di- 
vided into two groups, which we will call types A and 
B. The markings on the edges of the tiles indicate the 
matching rules - abutting edges of adjacent tiles have 
the same markings in a perfect tiling. It is readily seen 
that in an undefected tiling, a given row can consist of A 
or B type tiles only, with no mixing, and thus, we may 
characterize a row as being A-type or B-type. 

The KC tiling differs from other tilings studied in that 
it is not generated by recursive substitution (inflation) |9J , 
and by the fact that it has a ground state degeneracy 
with extensive entropy. This should be contrasted with 
tilings created from a simple inflation rule where the de- 
generacy scales as a power of the system size[ll] . All the 
ground states, however, are characterized by a quasiperi- 
odic arrangement of ^4-type and 5-type rows. The finite 
temperature behavior of this model is striking - as the 
temperature is increased from zero, the rows, still identi- 
fiable as A or B type, order periodically with decreasing 
period. At high enough temperatures, the rows lose their 
A or B identity, and the system becomes disordered. 
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The proof that there exists a zero-energy ground state 
is equivalent to showing that a perfect tiling exists. To do 
this [8 , we note that there are two numbers which char- 
acterize the n th row in a perfect tiling - the "frequency" 
a n and the number q n which are related through the 
mapping 



[_ma n J — L( m ~ 1 )r:t,i J 



ttfl+l = Qn OL r , 



(1) 



where 



l 
3 

1 < a r 



< 1 

< 2 



(2) 



Note that we shall take n to increase in the -y direction. 

The reason for dividing the tile set into the A and B 
groups becomes apparent if we denote the markings of 
the top, bottom, right, and left edges of a tile by {t,b,r,l}, 
as shown in Figure [2] When needed, we shall indicate 
the position of the tile as a subscript; thus t m n refers to 
the marking of the top edge of the tile centered at (m,n), 
etc. It is easily verified by inspecting Figure [T] that the 
markings on each tile satisfy the relation 



Qn ^m.n ~l~ ^r> 



(3) 



2 or |, re- 



with a tile of type A or B having q n 
spectively. The mapping a n +i = q n a n has no periodic 
points, since a n+q /a n — 2 p /3 q ~' p ^ 1 for any positive 
integers p and q; this implies that the resultant tilings 
are not periodic. The a n are distributed densely, but not 
uniformly, in the range ( | , 2) . 

To show that a perfect tiling exists, we must solve 
Equations [3] for all m,n with q n derived from Equa- 
tions [T] and [2] whilst demanding that b mn — t m n+ i and 
Im.n — ?*m-i,n- It is easily verified that the markings 

t m ,n = [ma n \ - [(m - l)a n J (4) 
r m ,n = Qn [ma n \ - [mq n a n \ 

constitutes a solution, where \x\ is the greatest integer 
less than or equal to x (see Figure [2| ; this is one example 
of a perfect tiling. 

To facilitate later discussion, we shall employ a useful 
mapping. Let us define the variable 4> n = nuj , where 
ujq = log2/log(6) w 0.3868. Decomposing <f) n into its 
integer and fractional parts gives 



<p n = [nu \ + {</>„} 



(5) 



where {4> n }, the fractional part of <j> n , is defined through 
its relation to a„ by 



log (a n ) + log (3) 
log (2) + log (3) 



(0) 




\mq rl Q n \ — l(m — Y)q n a., 



Figure 2: The marking scheme for the tile at (m,n). Note 
that n increases in the -y direction. 
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Figure 3: An example of two different 2x2 patches with the 
same outer markings; these may be exchanged with no match- 
ing rule violations. 



This maps a n € [|, 2) onto {4> n } € [0, 1), and allows us 
to obtain [L3] an explicit formula for q n : 
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Qn = 



-2 (LO + iHJ - L^oJ) + 2 



Sequences of this type are called Sturmian sequences, and 
are well known in the context of automatic sequences|3j. 
Here ujq is irrational, and this gives a quasiperiodic se- 
quence of the two "letters" 2 and | , with the consequence 
that in the ground state the rows appear in a quasiperi- 
odic sequence. 

As noted above, this ground state is only one of 
many, and in fact, there is an extensive ground state 
entropy [L3], This may be inferred from Figure [3] where 
two patches with the same markings on their exterior are 
presented (there exist larger patches with the same prop- 
erty as well). This means that starting from some given 
ground state, we may obtain another by randomly ex- 
changing the patches shown in Figure [3] (and any other 
pair of patches with the same exterior markings), pro- 
vided, as we have verified p"5]. that they appear with a 
finite density in the ground state. This implies that al- 
most all of the ground states are disordered in the sense 
that their patch entropy[12 scales as the patch size for 
large enough patches. This notwithstanding, in all the 
ground states, the rows are arranged in a quasiperiodic 
sequence of A-type and B-type, and it is also true that a n 
is unchanged for each row|15j. which is relevant to what 
follows. 
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Figure 4: Finite temperature configuration averaged over time 
exhibiting 3/8 periodicity. The row structure is evident, with 
colors as in Figure [I] Here L = 150, J x = J y , and T = 0.304. 




Figure 5: Qmax and Cy as a function of temperature, for a 
150 x 150 system with J x — J y . . 

Numerical studies of this model show that as the tem- 
perature is raised from 0, it goes through a series of phase 
transitions, where the A-B sequence of rows is periodic, 
and where the period decreases with increasing temper- 
ature. In analogy to the CI transition, we shall refer to 
these periodic phases as commensurate phases. In Fig- 
ure]^ we show a portion of a time averaged configuration 
from a 150 x 150 system (with J x = J y ) at T = 0.304 
(in units of J x ), which was obtained by parallel temper- 
ing, where the color coding is as in Figure [I] The A-B 
sequence is periodic with period 8, with 3 B rows per 
period. Characteristic defects are also present. At high 
enough temperature, of course, the rows lose their A or 
B character, and the system goes to a disordered phase. 

The different phase transitions are best traced by the 
winding number SI = AT N . B . T . Since at low T the rows 

° N A +N B 

are essentially pure type A or B, SI essentially counts the 
fraction of type B rows in the system. For an ensemble 
of systems at any given temperature, SI may take a va- 
riety of values, where the distribution has a well-defined 
maximum value, fl m ax, which is typically close but not 
identical to the average (S7) [21]. For large systems, both 



tlmax and (SI) will equal wo ~ 0.3868 at T = and ap- 
proach yjj ~ 0.46 as T — > oo. The first order nature 
of the transitions between the different commensurate 
phases can be verified by looking at the distribution of 
SI near the transitions, which exhibits two well separated 
peaks, with one overtaking the other as the transition 
temperature is crossed |15j . This suggests that 0, max (T) 
is a reliable indicator of these transitions. 

In Figure [5] we show curves of VL max and the specific 
heat Cy vs. T, for a 150 x 150 system. At low tempera- 
ture, the system is in its ground state, and S72 maa; = uiq (to 
within i). As temperature is increased, Sl max undergoes 
three jump discontinuities before becoming fully contin- 
uous as the rows are no longer homogeneous. While the 
step-wise behavior of £l m ax is a clear indication of the 
existence of the commensurate phases and the first or- 
der nature of the transitions between them, its value at 
the plateaus does not give the exact periodicities due to 
limited resolution and the existence of defects, and these 
then must be inferred by looking at the configurations 
themselves, as in Figure [4] The first order transitions 
are accompanied by small bumps in the specific heat at 
the same temperatures. In Figure [5] we identify phases 
with periods of 8 and 13 rows with the lower and middle 
plateaus, respectively. In our examinations of systems 
of sizes up to L — 300, we have identified phases with 
SI = |, |, jg, and |j, depending on system size and val- 
ues of the coupling constants. At these temperatures, 
we have verified numerically that the rows substantially 
maintain their pure A or B character. The broad peak 
in the specific heat at T ~ .37 attends the loss of purity 
in the rows. 

These results may be understood by an effective coarse- 
grained description of this system, appropriate for low 
temperatures, when the system may be considered as 
composed of A and B type rows|22 . This description 
resembles the Frenkel-Kontorova model[6 , and exhibits 
a competition between length scales. To see this, note 
that for the perfect tiling, the value of the frequency a n 
for the n th row may be calculated from the tile markings: 
a n = limj^oo ~ Z)ot=1 tmn - From tms i using Equations 
[5] and [6] we can compute <fi n . Now although defects enter 
the rows at finite T, we may still define two frequencies, 
using the markings of the top and bottom of a row. The 
"top frequency" is defined as — 4 Ylm=i ^ mn while the 
"bottom frequency" is given by — \ Sm=i ^mn- For 
a perfect tiling, — a^ +1 , and therefore, from Equa- 
tion [l] a^/a^ — q n , but this will typically not be the 
case for defected tilings. This notwithstanding, the <f> n 
at finite T may be inferred from using Equation [6] in 
the same manner as for the perfect tiling. These will be 
the variables used in our effective description. 

To construct an effective free energy for this model we 
assume that the dominant contributions come from the 
entropy of the rows, and the energy due to mismatches 
between the rows, each of which is characterized by its 
frequencies a T and a B . We argue that the energy cost 
associated with an imperfect interface between rows n 
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and n + 1 goes as L J y | — a^ +1 1 , where L is the length 
of the row. Clearly, this term is zero in the ground state, 
and numerical simulations bear out this functional form 
for low temperatures[15 . At this level of coarse graining, 
a row is characterized by its frequency a T (or equiva- 
lently a B ), and its entropy should be a function of this 
frequency which is extensive in L, so that we shall write 
the entropy of the n th row as Is \f*n) ■ 

Taken together, we get that the free energy is given by 

F / L = ELi J y - a n+i \ ~ T * («[)■ :t is conve- 

nient to express this in terms of the variables <j) n discussed 

above. The free energy is then of the form 

£ = ^2 O^n^m+l) - Ts (<j) m ) (7) 

m 

where g {<j) m , 4>m+i) is a function that favors (j) m +i = 
<p m + luq, which holds identically in the ground state. 
The entropy s (4> n ) depends only on the fractional part 
{4> n }i and thus it is a periodic function with period 
one. It is the competition between these two length 
scales which gives the novel behavior observed. Although 
it is tempting to expand g (</>„, 4> n +i) to first order as 
g (4> n , 4> n +i) oc \<p n +i — <i>n — <^o | , we note that such an 
expression fails when both {<j) n } and {<fi n+ i} are larger 
than 1 — ujq , since this would imply two adjacent B rows, 
which carries a disproportionately large energy cost. 

The equilibrium configuration of the <f> n can be ob- 
tained by minimizing F . As in the FK model, the first 
term favors an incommensurate phase with a winding 
number ujq while the entropy favors a commensurate 
configuration. The temperature T plays the role of the 
strength of the periodic potential, so that commensurate 
phases are expected at high temperature while incom- 
mensurate phases are expected at low temperatures. The 
commensurate phases that are expected are those with a 
winding number close to loq such as 1/3, 3/8, 5/13, etc. 
We have observed some of these phases in our numerical 
study, as seen in Figure [5] where their presence is indi- 
cated by the plateau values of Vl max . At still higher tem- 
peratures the segregation into type A and B rows breaks 
down, resulting in a disordered phase. It is interesting 
to speculate about the low T behavior of this system. It 
might be that only at T — an incommensurate phase 
appears, but it could be that such a phase, possibly with 
power-law correlations, is stable at finite T. These issues 
will be addressed in future work. 
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